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ON  THE  MICROWAVE  INTERACTION  WITH  MATTER 
AND  MICROWAVE  BREAKDOWN  OF  AIR 

1.  PROPAGATION  OF  ELECTROMAGNETIC  RADIATION  IN  AN  ATMOSPHERE 

1.1  Introduction 

The  propagation  of.  electromagnetic  radiation  in  the  atmosphere 
and  in  outer  space  is  of  considerable  interest  in  many  areas  of  applied 
physics.  The  most  prominent  research  field  for  radiation  propagation 
is  the  communication  sciences.  The  propagation  of  radiation  in  the 
atmosphere  can  be  quite  complex.  This  complexity  is  due  to  a  variety  of 
physical  processes  which  occur  during  the  interaction  of  radiation  with 
atmospheric  constituents.  Furthermore,  the  electromagnetic  radiation 
constitutes  a  large  spectrum  of  frequencies  which  range  from  x-rays  to 
radio  waves,  with  each  frequency  interacting  with  the  atmosphere  in  a 
different  manner. 

In  this  report  we  restrict  ourselves  to  microwaves  and  will  dis¬ 
cuss  some  of  their  applications  and  propagation  aspects  in  the 
atmosphere.  One  application  utilizes  microwave  radiation  to  interact 
with  a  target  in  the  atmosphere  or  in  outer  space.  Depending  on  the 
intensity  of  the  incident  radiation,  the  target  can  be  heated  to  a 
desired  temperature,  melted  or  vaporized.  These  results  can  be  accom¬ 
plished  by  focussing  the  microwave  radiation  on  the  target  of  interest. 
Furthermore,  microwave  radiation  can  be  focussed  on  a  given  area  in  the 
atmosphere  to  create  a  degree  of  ionization  which  depends  critically  on 
the  wavelength  of  the  incident  radiation. 
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The  creation  of  an  ionized  medium  in  the  atmosphere  has  several 
interesting  consequences.  One  of  these  is  that  the  ionized  air  acts  as 
a  reflector,  or  an  attenuator,  of  incoming  radiation.  Another,  is  that 
the  ionized  air  is  a  source  of  radiation  in  the  visible  and  infrared 
which  could  adversely  impact  sensors  operating  in  the  visible  or  the 
infrared. 

These  applications  depend  on  the  understanding  of  the  interaction 
of  the  microwave  radiation  with  targets.  This  interaction  determines 
the  radiation  intensity  required  for  effecting  a  desired  change.  This 
is  discussed  in  Chapter  2  of  this  report.  Transport  of  high  intensity 
radiation  requires  a  knowledge  of  the  air  breakdown  threshold  power. 

This  is  discussed  in  Chapter  3  of  this  report. 

1.2  Radiation  Intensity  and  Focussing 

The  applications  of  microwave  radiation  discussed  in  Section  1.1 
require,  in  general,  a  high  intensity  radiation.  One  of  the  well-known 
techniques  for  obtaining  high  intensity  radiation  is  to  focus  coherent 
radiation  on  a  small  area.  The  smallest  area  on  which  radiation  can  be 
focussed  is  determined  by  the  diffraction  effects.  Such  techniques  have 
resulted  in  laser  intensities  in  10  watts  per  cm  in  the  optical 
region.  Generally,  the  divergence  of  a  coherent  radiation  with  a  wave¬ 
length,  A,  emanating  from  an  aperture  of  diameter,  d,  is 

5  0*^  1,1 

If  the  focussing  distance  is  R^,  then  the  radius  of  the  spot  size  is 

rs  ~  7  Rf  1,2 

Thus,  for  example,  to  focus  the  microwave  radiation  over  an  area  with  a 

diameter  of  100  cm  at  a  distance  of  10  km,  Equation  1.2  requires  that 

the  antenna  aperture  be  10^  times  the  radiation  wavelength  (d  *  10^A) . 


The  approximate  relation  given  in  Equation  1.2  can  be  utilized  to 
calculate  the  desired  intensity  at  the  target  and  the  corresponding 
intensity  at  the  antenna.  In  a  more  realistic  maimer,  however,  the 
attenuation  of  the  radiation  in  the  atmosphere  should  be  considered 
along  with  the  intensity  distribution  on  the  focussed  spot. 

The  semi  quantitative  relation  of  Equation  1.2,  coupled  with  the 
approximate  relations  developed  in  Chapters  2  and  3  provide  a  coherent 
approach  to  the  understanding  of  certain  propagation  aspects  of  the 
microwave  radiation. 

1.3  Higher  Intensity  and  Coupling  Considerations 

The  interaction  of  the  microwave  radiation  with  a  target  and  its 
air  breakdown  are  two  aspects  of  interest  discussed  in  this  report. 

Both  problems  are  considered  on  theoretical  grounds  based  on  certain 
simplifying  assumptions.  However,  these  problems  require  experimental 
verification  as  well  as  more  detailed  and  self-consistent  theoretical 
and  computer  modelling.  These  requirements  are  necessitated  by  the 
lack  of  information,  especially  when  high  microwave  power  densities 
(P  >  106  watts)  are  considered.  Experimental  breakdown  studies  by 
pulsed  microwave  radiation  have  not  been  carried  out  for  any  power 
level  beyond  10  watts.  Higher  power  densities,  e.g.,  10  watts  are 
becoming  available  and  would,  in  principle,  enhance  the  local  effective 
electric  field  by  a  factor  of  30  which  clearly  will  change  the  free 
electron  distribution  and  hence  the  ionization  condition  and  the  emis¬ 
sion  characteristics  of  the  ionized  air.  In  the  area  of  the  interac¬ 
tion  with  a  target,  the  high  intensity  microwave  radiation  may  produce 
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a  vaporized,  ionized,  medium  on  the  target  surface.  Several  interest¬ 
ing  questions  arise  in  this  aspect.  For  example,  what  is  the  role 
of  this  atomized  region  in  terms  of  coupling  the  radiation  energy  into 
the  target?  Does  the  ionized  region  become  a  source  of  energy  that 
impacts  the  target?  To  answer  these  and  other  questions,  further 
experimental  and  theoretical  considerations  are  in  order. 


2.  IRRADIATION  OF  A  TARGET  BY  ELECTROMAGNETIC  RADIATION 


2.1  Introduction 

The  irradiation  of  a  target,  whether  a  conductor  or  a  dielectric 
material,  by  electromagnetic  radiation  will  raise  the  target  temperature 
and  may  crack,  melt  or  vaporize  it.  The  physical  changes  depend  on  the 
intrinsic  characteristics  of  the  material  and  the  amount  of  energy  it 
absorbs.  In  this  section  we  discuss  these  physical  changes  and  the 
amount  of  power  required  to  initiate  them.  This  can  be  accomplished 
by  calculating  the  target  temperature  as  a  function  of  time  and  its 
dependence  on  the  power  and  the  wavelength  of  the  incident  radiation. 

The  target  temperature  is  obtained  from  an  analytic  solution  of  the 
heat  equation  in  one  dimension.  From  this  analytic  solution,  approxi¬ 
mate  expressions  for  the  surface  temperature  are  obtained  and  are 
utilized  for  quantitative  estimates  of  the  power  required  to  vaporize 
a  target. 

When  the  radiation  is  incident  on  a  target,  a  certain  amount  of 
the  energy  is  reflected.  The  reflectivity  of  the  target,  obviously, 
influences  the  total  incident  radiation  required  to  vaporize  a  target. 

In  this  chapter  we  give  a  simple  relation  for  the  reflectivity  of  good 
conductors  as  a  function  of  the  wavelength  of  the  incident  radiation. 
This  will  be  incorporated  in  the  approximate  target  temperature 
equation  to  give  the  total  power  necessary  to  induce  the  physical 
changes  on  a  target. 

2.2  Heat  Equation  and  the  Target  Temperature 

The  equation  describing  the  heat  flow  into  a  solid  target  is1 

T7  r?2  t  j.  n  -  _  r  ^T  /■- 


where  K  is  the  thermal  conductivity  of  the  target,  p  its  density,  Cp  is 
its  specific  heat  at  constant  pressure,  Q  is  the  rate  of  the  energy  in¬ 
put  per  unit  volume  and  T  is  the  temperature.  The  above  equation  can  be 
rewritten  as 


v  2T  + 


Q,  1  dT 
T  ~  k  at 


(2-2) 


where  k  is  the  thermal  diffusivity  of  the  target  and  is  defined  as 
which  is  the  ratio  of  the  thermal  conductivity  to  the  heat  capacity  of 
the  target. 

Consider  the  solution  to  Equation  (2-2)  for  a  uniform  irradia¬ 
tion  incident  perpendicular  to  the  surface  of  a  semi- infinite  slab  of 
material.  We  seek  a  solution  in  one  dimension  for  an  incident  power 
which  is  absorbed  as  it  traverses  the  target  according  to  Q=Q0  e^*X 
where  a  is  the  absorption  coefficient  of  the  target  whose  surface 
boundary  is  at  x=o.  Assuming  that  the  target  temperature  is  T*o  at  t=o, 
the  solution  to  Equation  (2-2)  is 


T(t,x)  *  ol 


+  2o— 
2K  o? 

+  3cl_ 
2K  o? 


e-« 


o?  kt-otx 


erfc  £a(kt) 


2(H)' 


b 


e<#kt«x  erfc  |a(kt)>,  . 


where  Qq  is  the  absorbed  power  and 
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i  erfc  a  =  J"  erfc  y  dy 


00  O 

2  r»  ~  y 

erfc  a  =  —  J  e  dy 
a 


i  erfc  o  =  7=- 
i  erfc  co  =  o 


The  surface  temperature  of  the  target  can  be  obtained  immediate¬ 
ly  from  Equation  (2-3)  by  setting  x=o,  that  is; 


Ct.ol  ■  1  .  f*.  eric 


(2-4) 


The  solution  to  Equation  (2-2)  was  obtained  under  the  assump¬ 
tion  that  the  thermal  conductivity  and  the  theimal  diffusivity  of  the 
target  are  constant  with  temperature.  This  is  reasonable  for  metals 
where  these  coefficients  change  slightly  with  the  target  temperature. 
Table  1  gives  the  thermal  conductivity  and  the  heat  capacity  of  copper 


as  a  function  of  temperature  . 


Table  1 


K  /Watt  \ 
ycm  °c  J 

K/k  /  Joule  \ 
l  mole  0 c  J 


0  °c 

100  °c 

300  °c 

u 

o 

o 

o 

3.85 

3.82 

3.76 

3.5 

24.6 

24.9 

25.2 

27.7 

However,  the  thermal  conductivity  is  reduced  by  a  factor  of  two  when 


the  solid  is  melted. 

2.2  The  Reflectivity  of  the  Target 

In  Equations  (2-3)  and  (2-4),  Q0  represented  the  power  absorbed 
in  the  target  which  generates  the  heating  and  the  rise  in  the  target 
temperature.  However,  when  radiation  is  incident  on  any  surface,  part 
of  the  radiation  is  transmitted  while  another  part  is  reflected.  There¬ 
fore,  it  is  of  interest  to  know  the  reflectivity  or  the  transmittance  of 
the  target  for  each  wavelength  in  order  to  obtain  the  total  incident 
radiation  on  the  target.  Therefore,  if  we  defined  P  as  the  total  inci¬ 
dent  power,  then 

Qo  =  U-R)  P  (2-5) 

where  R  is  the  reflectivity  of  the  target. 

Consider  the  case  of  normal  incidence  of  electromagnetic  radia¬ 
tion  on  the  surface  of  a  good  conductor.  The  transmittance  of  a  good 
conductor  is  related  to  the  skin  depth,  5,  asJ 

1_R  *  "T  6  (2-6) 

where  A  is  the  wavelength  of  the  incident  radiation  and 

®  =  /ZniJic  cr  (2-7) 

In  Equation  (2-7),  n  is  the  permeability  and  a  is  the  electrical  con¬ 
ductivity  and  u :  is  the  angular  frequency  of  the  incident  radiation. 


Using  Equation  (2-7)  into  (2-6),  one  obtains 


1-R  =  3.46  X  loP  -  (2-8) 

//i  A  Q 

In  general,  the  electrical  conductivity  of  a  good  conductor 4 
is  cr  ~  1017  sec"1  while  the  permeability  is  ~  1  (for  non  magnetic  mate¬ 
rial).  For  radiation  at  A=3  cm,  the  amount  absorbed  in  a  conductor  is 
~  6  X  10"4  .  In  Figure  2.1  are  shown  1-R  values  for  copper,  Cu,  and 
aluminum,  Al,  over  a  wavelength  range  0.01-3  cm.  For  example,  at 
A=3  cm,  the  amount  of  radiation  transmitted  into  Cu  and  Al  are 
2.3  X  10"4  and  3.4  X  10  11  ,  respectively.  On  the  other  hand,  for  car¬ 
bon,  the  transmittance  at  A=3  cm,  is  8  X  10~3  which  is  larger  by  an  or¬ 
der  of  magnitude  compared  to  that  for  Cu  and  Al.  This  is  because  the 
electrical  conductivity  of  carbon4  is  6.3  X  1014  sec’1  which  is  lower 
by  two  orders  of  magnitude  compared  to  good  conductors. 

For  dielectric  materials  on  the  other  hand,  assuming  normal  in- 

c 

cidence,  the  reflectivity  is 

r=(htt)2 

where  n  is  the  index  of  refraction.  Dielectrics  such  as  epoxy  cast 
resin,  polyethylene  and  rubber,  have  index  of  refractions6  of  n=1.7  for 
the  visible  light.  For  longer  wavelength  (A=10  cm),  their  refractive 

*7 

indices'  range  from  1.4  to  1.7  implies  small  reflectivities.  This 
implies  that  in  a  sense,  less  total  incident  power  will  be  required  to 
crack,  heat  or  vaporize  dielectric  materials. 
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2.3 


Vaporization  of  a  Solid 

In  Section  2.2  an  analytic  solution  was  given  for  the  tempera¬ 


ture  at  any  point  inside  a  solid  and  as  a  function  of  time  (see  Equa¬ 
tions  2-3  and  2-4).  These  equations  can  be  utilized  to  heat  a  target  to 
any  desired  temperature  knowing  the  power  of  the  incident  radiation. 
However,  certain  approximate  expressions  for  the  target  temperature  can 
be  obtained.  From  these  expressions  one  can  obtain  quantitative  esti¬ 
mates  for  the  power  needed  to  melt  or  vaporize  a  target. 

Consider  the  surface  temperature  of  a  target  as  given  by  Equa¬ 
tion  (2-4).  An  approximate  expression  can  be  obtained  for  early  times 
when  o?  kt  «  1  or  o?  k  «  ^-.  In  this  case , 

T  (t,o)  -  (1-R)  o.  |P.t,  (2-10) 

where  use  is  made  of  Equation  (2-5).  The  concept  of  an  early  time  can 
be  clarified  if  one  defines 

to  *  SHc  “  F"  (2~n) 

where  the  absorption  length  is  of  the  order  of  the  skin  depth  and  that 
tg  is  the  time  required  for  heat  to  diffuse  a  distance  of  <$.  Thus 
for  any  time  t  <  t^,  Equation  (2-10)  can  be  used  to  calculate  the 
power  needed  to  raise  the  target  surface  to  melting  or  vaporization 
temperatures.  The  skin  depth  for  good  conductors  is  vio'^cm  at  \  -  5 
cm  and  k  1  cm2/ sec,  thus  t^  -  10  ^sec.  However,  of  more  interest 
is  heating  a  target  for  longer  times. 
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Thus  we  would  like  to  obtain  an  approximate  relation  from  (2-3)  and 
(2-4)  for  times  t  where 


tc  <  t  <  t 


u 


(2-12) 


,2 

Here  t,  =  ^  with  x  being  the  target  thickness 
large  compared  to  the  skin  depth  (2,>>6).  The 
here,  implies  that  heat  has  diffused  beyond  a 
temperature  of  the  back  surface  of  the  target 
approximation  yields 


which  is  usually  very 
concept  of  a  long  time, 
skin  depth  but  that  the 
is  unchanged.  Such  an 


T(0,t) 


2(l-R)P^kt^ 


(2-13) 


which  we  will  use  to  estimate  the  total  power  needed  to  vaporize  a 
skin  depth  of  the  target  material.  This  is  done  by  ignoring  the  liquid 
phase  of  the  solid  which  is  valid  since  the  latent  heat  of  fusion  is 
much  smaller  than  the  latent  heat  of  vaporization.  The  heat  of  fusion 
is  also  much  smaller  than  the  quantity  of  heat  needed  to  raise  the  tar¬ 
get  to  vaporization  temperature. 

The  time  to  reach  vaporization  temperature  for  any  power  level 
can  be  obtained  by  rewriting  Equation  (2-13),  that  is 


_  m  K  2  Tv2  l 

v  "  4  k  (1-R)2  F* 


(2-14) 


It  is  obvious  from  Equation  (2-14)  that  higher  power  levels  are  needed 
if  one  desires  shorter  times  to  reach  vaporization  temperatures. 


2.3.1  Applications  to  Good  Conductors 

Consider  copper  where  K  =  3.85  watt  cm"1  °c  ,  k  -  1.1  cnf  sec*1 
and  Ty  is  2300  °c6.  The  time  to  reach  vaporization  temperature  is 

*  _  5.5  X  107 
v  ~  (1-R)2  ? 

where  P  is  in  watts  cm"2  .  This  relation  is  illustrated  in  Figure  2.2 
for  several  wavelengths  of  the  incident  radiation.  It  is  clear  from 
this  Figure  that  vaporization  temperature  is  reached  faster  and  at  a 
lower  power  level  for  shorter  wavelengths  of  the  incident  radiation. 

Another  interesting  application  of  Equation  (2-13)  is  when  the 
incident  power  duration  is  specified  and  one  is  then  interested  in  the 
power  needed  to  vaporize  the  solid.  Using  the  values  of  K,  k  and  Tv 
given  earlier  for  copper.  Equation  (2-13)  reduces  to 

P(l-R)  =  7'-5- *  ^  (2-15) 

/t 

Now  consider  an  incident  radiation  at  A=3  cm  and  t=10"2  sec,  with  1-R 
being  2.7  X  10"“*  one  needs  P=2.8  X  l(f  watts/cnf  to  heat  copper  to  its 
vaporization  temperature.  This  is  equivalent  to  2.8  X  1012  watts/nf  . 

The  analysis  leads  to  P=1.5  X  1012  watts/n?  for  A1  where  K  *  2.37  watt 
cm*1  °c  ~l ,  ks1.0  ai?  sec"1  and  Tv=2460  °c^«  These  power  levels  are  in 
agreement  with  previous  calculations  reviewed  in  Reference  (7) . 

Once  the  vaporization  temperature  is  reached,  the  rate  of  the 
material  removal  and  the  vaporized  depth  of  the  material  can  be  esti¬ 
mated.  The  rate  of  vaporization  or  the  surface  ablation  velocity,  V  , 

”  cl 
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can  be  expressed  ,  as 


ya  .  P 


where  heat  losses  due  to  convection  and  conduction  are  ignored.  In 
Equation  (2-16),  p  is  the  density  of  the  solid,  Ly  is  the  latent  heat 
of  vaporization.  Thus  the  vaporized  depth,  i  ,  is  obtained  from 
Equations  (2  -  16  )  for  the  velocity  and  Equation  (2-14)  for  the 
vaporization  time, 

V  =  Vv  C2'17) 

For  copper,  Ly  =  4000  J/gm,  p  =  7  gm/cm3  and  Cp  -  0.09  cal  gm'1  c°*. 

Using  these  values  into  Equations  (2-16)  and  (2-17)  with  a  power  of 

2.8  x  10  watt/ cm  at  X  =  3  cm,  one  obtains  i.  »  0.025  cm.  Thus  to 

vaporize  one  centimeter  depth  of  a  copper  target,  one  needs  40  pulses 
-2 

of  10  sec  duration  for  each  pulse.  This  estimate  assumes  that  the 
energy  is  coupled  to  the  target  in  the  same  way  for  each  subsequent 
pulse.  Such  an  estimate  may  be  appropriate  for  a  moving  target.  For 
a  still  target,  however,  the  amount  of  surface  material  vaporized  after 
the  first  pulse  may  act  as  an  absorber  of  the  incident  energy  in  the 
pulses  that  follow.  This  may  decouple  the  incident  radiation  from  the 
target.  However,  further  heating  of  the  vaporized  region  could  result 
in  a  plasma  whose  radiation  in  shorter  wavelengths  is  absorbed  by  the 
target.  Obviously  the  problem  is  rather  complex  and  requires  a 
comprehensive  treatment. 

2.3.2.  Application  to  Carbon 

14  -1 

The  electrical  conductivity  of  carbon  is  6.3  x  10  sec 
which  is  lower  by  two  orders  of  magnitude  compared  to  that  for  good 
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metallic  conductors.  This  characteristic  makes  carbon  a  better  absorb¬ 
er  of  radiation  compared  to  good  conductors.  For  example,  1-R  for  car¬ 
bon  is  8  X  10~3  at  A=3  cm  compared  to  2.7  X  10"4  for  copper.  For  this 
reason  it  is  of  interest  to  apply  the  approximate  vaporization  relations 
of  Section  2.3  to  carbon. 

Using  Equation  (2-13)  for  carbon  where  K  =  0.24  watt  cm'1  °c  ~l  , 
k-0.16  cn?  sec’1  and  Tv  =  4830  °c  we  obtain 


p  *  2-6  X  1C? 

(1-R  Vt 


(2-18) 


which  yields  P=3.2  X  106  watts/cn?  for  X=3  cm  and  1-R=8  X  10  3  .  This 
value  of  power  to  raise  carbon  to  its  vaporization  temperature  is  small 
by  two  orders  of  magnitude  compared  to  that  for  copper. 


2.3.3 


jlication  to  Dielectric  Material 


The  characteristics  of  the  dielectric  materials  are  different 


than  those  of  conductors.  Therefore,  it  is  of  interest  to  consider 
these  materials  under  the  impact  of  electromagnetic  radiation.  To  apply 
the  equations  developed  earlier  for  solids,  one  needs  the  amount  of  ra¬ 
diation  transmitted  into  the  solid,  i.e.  the  quantity  1-R.  For  the  di¬ 
electric  materials  the  reflectivity  for  normal  incidence  is  defined  by 
Equation  (2-9)  which  depends  on  the  index  of  refraction.  In  Table  2, 

g 

the  dielectric  constant  for  several  plastics  are  given  as  a  function 
of  the  radiation  frequency. 
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Table  2 


^(Hz) 


Plastics 

103 

106 

lcf 

3X1  CP 

2.5X10 

Epoxy  Resin 

3.67 

3.62 

3.35 

3.09 

Polyethylene 

2.26 

2.26 

2.26 

2.26 

2.26 

Polyhexamethylene 
-  adipamide  (Nylon) 

3.5 

3.14 

3.0 

2.84 

2.73 

Polytetrafluoroethylene 

(Teflon) 

2.1 

2.1 

2.1 

2.1 

2.1 

For  the  region  of  interest,  the  dielectric  constants  for  these 
plastics  change  very  slightly  with  the  radiation  frequency.  Therefore, 
the  index  of  refraction  is  also  constant  with  respect  to  the  incident 
wavelength.  The  corresponding  values  of  n  for  polyethylene  and  teflon 
are  1.5  and  1.45,  respectively.  This  implies  that  1-R  for  these  two 
plastics  are  0.96  and  0.97,  which  means  that  they  are  practically  trans¬ 
parent  to  microwave  and  radio  frequencies.  Table  3  gives  the  necessary 

9 

data  on  polyethylene  and  teflon. 


Table  3 


Density 

Thermal  Conductivity 

Specific  Heat 

gm/cm3 

10"4  cal/sec/ar?  /l°c  /cm 

cal/° c  /gm 

Polyethylene 

1.28 

3.8 

0.32 

-  oxide  based 

Teflon 

2.17 

6.0 

0.25 

These  data  can  be  used  to  obtain  the  thermal  diffusivities 
which  are  9.2  X  10"4  aif  /sec  and  1.1  X  10'3  aif/sec  for  polyethylene 
and  teflon,  respectively.  With  the  information  in  Tables  2  and  3  on  the 


transmittance,  thermal  conductance  and  thermal  diffusivity,  we  can  esti 
mate  the  power  needed  to  decompose  these  dielectrics.  However,  one 
must  estimate  the  absorption  coefficient  or  the  skin  depth  in  order  to 
apply  the  approximate  relations  given  in  Equations  (2-10)  and  (2-13). 
The  absorption  coefficient  for  teflon  at  A  *  0.005  cm  is10  200  cm"1 
which  implies  a  skin  depth  of  5  x  10 an.  The  time  for  heat  to  dif¬ 
fuse  a  skin  depth  is  0.025  sec.  Thus  one  can  utilize  Equation  (2-13) 
for  times  longer  than  0.025  sec,  say  2  x  10  ^sec.  Using  the  relevant 
data  for  teflon  into  Equation  (2-13),  where  the  decomposition  tempera¬ 
ture  is  35 08,  one  obtains 


This  implies  that  for  a  pulse  duration  of  10_1sec  one  needs  =47  watts/ 

cm"  to  decompose  teflon  with  a  radiation  at  A  =  5  x  10'^cm.  At  longer 

wavelengths,  however,  the  absorption  coefficients  are  smaller  and  can 

be  obtained  if  the  loss  angle  or  power  factor  are  known.  The  absorp- 

12 

tion  coefficient  can  be  expressed  as 

a  a  ^  n,  tan  5  *  0.0001  to  0.05 
and 

2tt  "  ""  ,,vm' 

a  3  ~xn  vl  ♦  tan  o  -  1  ,  tan  5  *  0.05  to  50 
2 

where  tan  6  is  the  loss  factor.  For  polyethylene  the  loss  factor  is  * 
6  x  10  for  A.  from  3  an  to  10  cm.  This  implies  that  the  absorption 
coefficient  is  small  and  the  skin  depth  large  in  these  wavelength 
ranges  and  one  must  select  the  appropriate  time  and  Equations  to 
estimate  the  power  needed  to  decompose  these  dielectrical  material. 


2.3.4.  Remarks  on  Target  Heating 

When  a  target  surface  is  irradiated  by  a  sufficient  flux  from 
a  radiation  source,  certain  amount  of  the  surface  material  is  vaporized. 
This  vaporized  layer  has  an  atomic  nature  and  its  absorbing  character¬ 
istics  are  very  different  than  those  of  a  solid  and  a  polished  surface. 
The  vaporized  layer  absorbs  the  incoming  radiation  and  may  decouple  it 
from  the  target  surface.  Thus  the  calculations  presented  here  are,  in 
principle,  applicable  to  a  moving  surface  whereupon  the  vaporized  layer 
does  not  pose  a  decoupling  problem.  However,  a  comprehensive  treatment 
requires  a  detailed  time  dependent  model  calculation  which  takes  into 
consideration  the  vaporized  layer  and  its  effects  on  the  heating  and 
vaporization  of  target  material  as  well  as  the  hydrodynamic  effects. 


le  transmittance  of  aluminum  and  copper  as  a  function  of  wavelength 


3.  AIR  BREAKDOWN 


3.  Introduction 

In  Section  Two  of  this  report,  quantitative  relations  and  esti¬ 
mates  were  given  for  heating  and  vaporization  of  solid  materials  under 
the  impact  of  electromagnetic  radiation.  However,  radiation,  generally, 
propagates  in  a  medium,  such  as  air,  before  reaching  the  point  of  inter¬ 
est.  Propagation,  especially  in  air,  is  accompanied  by  several  physical 
processes  which  affect  adversely  the  transportation  of  radiation.  One 
of  the  important  processes  that  limits  the  propagation  of  radiation  in 
air  is  the  air  breakdown  which  will  be  discussed  in  this  section.  This 
discussion  will  provide  quantitative  relations  which  determine  the 
breakdown  threshold  as  a  function  of  the  power  and  the  wavelength  of 
radiation.  Such  information  is  essential  in  order  to  avoid  air  break¬ 
down  when  radiation  is  transported  and  focussed  at  any  desired  point  in 
space. 

3.2  Air  Breakdown  by  Electromagnetic  Radiation 
3.2.1  Theory  of  Breakdown 

The  breakdown  in  air  or  any  other  gaseous  medium  is  exemplified 
by  the  creation  of  a  plasma.  This  occurs  when  sufficient  power  from  an 
electromagnetic  source  is  focussed  in  air,  whereupon  the  free  electrons 
present  gain  energy  under  the  influence  of  the  intense  electric  fields 
associated  with  the  focussed  radiation.  The  energy  gained  by  electrons 
is  expended  in  producing  more  electrons  by  electron  impact  ionization 
culminating  in  cascade  ionization  of  the  gaseous  medium. 

There  have  been  considerable  studies  of  breakdown  in  many  gase¬ 
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ous  elements  and  in  air  under  the  influence  of  microwave  radiation^' l3"?-5 
With  the  advent  of  lasers  and  high  power  Q- switched  lasers,  the  air  and 
gaseous  breakdown  studies  at  optical  frequencies  have  increased  consid¬ 
erably16,17. 

The  mechanism  for  gas  breakdown  can  be  explained  by  the  classi- 
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cal  microwave  breakdown  theory  ,  where  the  electron  gains  energy  from 

the  electric  field  only  when  it  collides  with  an  atom  and  the  sudden 

change  in  its  velocity  results  in  the  transfer  of  the  oscillation  energy 

into  random  motion.  This  energy  gain  process  has  been  shown  to  corres- 
18 

pond  to  the  quantum  mechanical  description  of  the  energy  gain  by  the 
free  electrons  from  the  radiation  field  through  the  free-free  transi¬ 
tions  in  the  field  of  an  atom. 

The  rate  of  energy  gain  by  the  free  electron  under  the  influ¬ 
ence  of  the  electric  field,  E,  of  the  incident  radiation  is 


df  \  _  e2  E2  v  coz 

mu)2  o;2+  v 


(3-D 


where  w  is  the  natural  frequency  of  the  oscillatory  electric  field,  vQ 

is  the  collision  frequency  of  the  electron  with  the  neutral  atoms.  On 
the  other  hand,  a  simplified  model  for  the  rate  of  energy  loss  in  air  is 


-  D  T  -  4B 
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Where  D  is  the  diffusion  coefficient,  A  the  diffusion  length,  e  is  the 
average  electron  energy,  v.  ,  v  and  v  are  the  collision  frequencies  for 
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ionization,  attachment  and  excitations,  respectively,  while  e ^  and  ee 
are  the  corresponding  energies  for  ionization  and  excitations,  respec¬ 
tively.  The  second  term  on  the  right  hand  side  of  Equation  (3-2)  re¬ 
presents  the  energy  loss  by  electrons  through  elastic  collisions  with 
the  neutrals.  Therefore,  the  net  energy  gained  by  electrons  is 

Furthermore,  the  time  development  of  the  electron  density  is  described 
by 


dne  D 

~3t  =  V'i  '  ne^a  '  O4) 

where  it  is  obvious  that  electron  attachment  and  diffusion  will  control 
th*»  cascade  ionization  process  and  must  be  overcome  for  breakdown  to 
occur.  Two  limiting  cases  of  Equation  (3-4)  can  be  of  great  interest. 
These  are  the  diffusion  controlled  and  attachment  controlled  cases  of 
air  breakdown.  At  atmospheric  pressure,  however,  the  air  breakdown  is 
attachment  controlled. 

The  condition  for  a  cw  breakdown  is  =  0  which  implies  that 


+ 


D 


The  diffusion  coefficient  can  be  expressed  as 
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(3-5) 
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Dp  *  3.2  X  lcfe 


(3-6) 


where  p  is  the  pressure  in  mn  Hg  and  e  is  the  average  electron  energy  in 

units  of  eV.  Thus  for  typical  energies  of  few  eV  and  atmospheric 

2 

pressures  of  760  iron,  the  diffusion  coefficient  is  D  “  8  X  l(f  (sec-*  ’ 

This  implies  that  the  rate  of  diffusion  from  an  area  of  1  cm2  is 
~  103  sec"1  which  is  very  small  compared  to  the  attachment  rate  of  elec¬ 
trons  to  0  .  This  latter  rate  is  ~  3  X  107  sec'1.  However,  at  higher 
2 

altitudes  or  lower  air  pressures,  diffusion  becomes  important  and  in 
certain  cases  even  a  dominant  controlling  factor  in  air  breakdown. 

3.2.2  Simple  Air  Breakdown  Calculations 

The  air  breakdown  calculation,  in  principle,  can  be  carried  out 
by  solving  Equation  (3-3)  and  (3-4)  with  the  knowledge  of  the  electron 
velocity  distribution  and  the  appropriate  atomic  processes  whose  rates 
are  electron  energy  dependent.  However,  reasonable  estimates  (lower 
limits)  can  be  obtained  for  the  breakdown  threshold  power  using  approx¬ 
imate  solutions.  The  rate  of  energy  growth  for  the  electron  in  the  ra¬ 
diation  field,  Equation  (3-1),  can  be  written  in  terms  of  the  incident 
radiation  flux,  I,  as 


d£  _  4ff#I 
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where  ^  ^  was  utilized  (E  =  27/1  (V/an)  with  I  in  watts/cnF  .  The 


threshold  condition  is  obtained  using  Equation  (3-7)  assuming  that  en¬ 
ergy  gained  by  the  electron  is  expended  to  ionizations  and  that  other 
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loss  mechanisms  are  ignored.  Thus  the  energy  gained  by  the  electron 
during  the  radiation  pulse,  At,  is  equal  to  the  ionization  potential  of 
the  neutral,  1^,  times  the  critical  number  of  generations,  K^, 


(3-8) 


Thus  the  threshold  power  for  air  breakdown  is 
I  •  K  (to2  +  v  2  ) 

1  CT  ^  C  J 

Tth  -  ■■oTirst— —  (3*9) 

For  a  pulsed  breakdown  the  first  electron  generates  another  and  then 
these  generate  two  more  electrons  and  so  on.  The  growth  of  the  elec¬ 
tron  density  is  described  by 


Ne  =  Ne(0)  eyAt  =  Ne(0)  2Kcr 


(3-10) 


Where  v  is  the  net  ionization  rate  and  Ng(0)  is  the  initial  density  of 

the  free  electrons,  generally  ~  1  cm"3.  Furthermore,  it  is  assumed 

that  at  the  end  of  the  radiation  pulse  the  electron  density  reaches  its 
cr 

critical  value,  N  ,  where 
e 


K 


cr 


=1.45  loge 


(3-11) 


The  critical  electron  density  depends  on  the  radiation  wavelength  that 
causes  the  air  breakdown  and  is  given  by 


dfe-n.ijr 
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(3-12) 


e  A2 

which  is  shown  on  Figure  3.1  for  several  wavelengths. 

To  estimate  the  threshold  power  for  air  breakdown  under  optical 
frequencies,  Equation  (3-9)  can  be  utilized  with  Kcr  values  ranging 
from  36  to  40.  Using  K  =  40  and  1^  =  14.8  eV  for  the  average  ioniza¬ 
tion  potential  of  air  molecules  in  Equation  (3-9)  results  in 


lth  =  8,6  X  10  16  (watts/cm2  )  (3-13) 

c 
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The  collision  frequency  can  be  estimated  to  be  v  *  5.3  X  lCr  p  where 
p  is  the  air  pressure  in  Torrs,  therefore,  Equation  (3-13)  reduces  to 


8.6  X  10~>6 
At 


o.2  +  (5.3  X'  1C?  p)2 
5.3  X  Ilf  p 


(3-14) 


For  optical  region  1016  which  is  much  larger  than  v  ,  one  obtains 

for  air  at  atmospheric  pressure  a  breakdown  threshold  power  of  2  X  10^ 

(watts/cm2  )  using  a  pulse  with  duration  of  10~6  sec.  This  value  is  of 

the  same  order  but  lower  by  a  factor  of  three  compared  to  the  breakdown 
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threshold  power  obtained  experimentally  ,  indicating  that  the  classi¬ 
cal  approach  can  explain  the  breakdown  process.  Thus  the  breakdown 
calculations  require  the  interrelation  of  many  factors  which  are  the 
ionization  rate,  the  energy  loss  rate,  the  average  electron  energy,  the 
electric  field  and  finally,  the  threshold  power  for  breakdown.  This 
requires  a  more  detailed  analysis  which  is  carried  out  in  the  next 
section  for  microwave  breakdown  of  air. 
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3.2.3  Calculation  of  Air  Breakdown  by  Microwave  Radiation 
21  22 

Approaches  '  to  air  breakdown  calculations  can  start  with  the 
knowledge  of  the  ionization  rate  of  air,  which  generally  is  ex¬ 
pressed  as  !ii.  .  This  quantity  depends  on  the  electron  energy  which  in 
P  £ 

turn  depends  on  the  critical  electric  field,  —  .  For  a  Maxwellian 

p 

electron  velocity  distribution  and  for  an  average  electron  energy  of 


e  ^  20  eV,  the  following  relation  for  —  has  been  obtained 


22 


^  -  2.5  X  107  [5.6  (O*5  +  0.6  (€)Hj  Exp  (3-15) 


where  it  is  assumed  that  in  air  ionization  is  basically  from  oxygen  mo- 
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lecule  (Ii  =>-12.5  eV) .  Furthermore,  an  empirical  relation  has  also 

Ee 

been  found  between  e  and  —  where 

P 


e  -  0.0675 


(3-16) 


This  relation  utilized  into  (3-15)  yields 


v . 

—■  =  2.5  X  107 


(3-17) 


which  is  shown  in  Figure  3.2.  Thus  the  problem  is  a  simple  one  for  air 
breakdown  when  it  is  attachment  controlled.  For  cw  the  condition  for 


breakdown  is 


or  equivalently 


v. 

i 

P 


(3-18b) 


a. 

The  value  of  —  is  known  and  thus  one  can  obtain  from  Equation  (3-17). 

However,  for  pulsed  breakdown,  the  ionization  rate  must  be  high  enough 
so  that  the  electron  density  reaches  its  critical  value  within  the  dura¬ 
tion  of  the  radiation  pulse.  This  requires  that 


v.  v  ,  N  cr 

i  .  a  _  1  .  e 

p  p  “  pAt  ge  Ne  CO) 


or 
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where  we  have  utilized  Equations  (3-10  and  (3-12).  Thus  for  any  wave- 

"i 

length,  At  and  p  we  can  obtain  —  •  From  Figure  3.2  or  Equation  (3-17) 

E 

we  obtain  the  corresponding  ~  and  hence  the  threshold  breakdown  power. 

Using  X  =  3  cm,  p  =  760  nri,  At  =  10'6  and  N  (0)  =  1  in  Equation  (3- 19b) 
Ee 

we  obtain  —  -  30  volts/cm  Torr. 

An  application  of  Equation  (3- 19b)  with  Ne(0)  *  1  for  a  fixed 

pulse  duration  (At=10"6sec)  is  shown  in  Figure  3.3  for  several  differ- 

v. 

ent  wavelengths.  In  this  Figure  -j—'s  are  given  as  a  function  of  air 


pressure  (altitude).  The  corresponding  effective  electric  fields  as  a 
function  of  pressure  (altitude)  are  shown  in  Figure  3.4.  In  this  figure 
the  independence  of  the  electric  field  on  the  radiation  wavelength  with 
increasing  pressure,  is  a  manifest  of  the  attachment  controlled  domain 
for  breakdown. 

3.2.4  Air  Breakdown  Calculation  With  Diffusion  Effect 

The  calculation  in  Section  3.2.3  of  this  report  was  carried  out 
ignoring  the  diffusion  effect.  However,  when  the  density  of  air  is  re¬ 
duced,  diffusion  may  become  an  important  factor  in  air  breakdown  and, 
therefore,  it  has  to  be  considered.  In  this  section  we  present  air 
breakdown  calculations  by  including  the  diffusion  effects  phenomenologi¬ 
cally.  In  addition,  we  present  the  effect  of  the  radiation  pulse  width 
on  the  breakdown  thresholds. 

The  criterion  for  a  cw  breakdown  is 


which  is  equivalent  to  Equation  (3-5).  For  a  pulsed  condition,  however, 
the  criterion  for  breakdown  is 

F ("i  -  "a  -  tf)  ’  ESt  loSe(NeCr/Ne(0)) 

This  implies  that  two  other  parameters  have  to  be  known  in  order  to 
calculate  the  critical  field.  These  are,  D  and  A.  The  diffusion  co¬ 
efficient  is  related  to  Ee/P  according  to 
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DP  =  29.0  +  0.9  j  X  10H 


(3-22) 


and  A  is  generally  a  characteristic  length  of  the  breakdown  geometry  and 
can  be  estimated.  In  Figure  3.S  we  present  breakdown  threshold  power 
calculations  as  a  function  of  pressure  for  several  wavelengths  of  the 
incident  radiation.  These  calculations  were  carried  out  with  A=0.75  cm 
and  At=  3X10 '6  sec.  It  is  obvious  from  this  Figure  that  the  threshold 
breakdown  power  increases  with  increased  frequency  of  the  incident  radi¬ 
ation.  However,  these  calculations  were  repeated  for  a  shorter  pulse 
width  i.e.  At=3X10~7  sec  to  demonstrate  the  effect  of  shorter  pulse 
duration.  These  latter  results  are  shown  in  Figure  3.6  where  it  is 
obvious  that  the  reduction  in  the  radiation  pulse  width,  for  each  wave¬ 
length,  results  in  higher  breakdown  threshold  power. 

A  final  remark  is  in  order  in  regards  to  breakdown  calculations. 

Free  electrons  in  ionized  air  decay  through  recombination  in  addition 

to  attachment  and  diffusion  from  the  ionized  region.  For  ionization 

by  microwave  radiation  (A  »  1  cm) ,  electron  densities  are  ng  <  1013 

cm'3  where  recombination  will  proceed  through  the  process  of  dissocia- 

22 

tive  recombination  which  depends  on  the  electron  temperature  .  The 
dissociative  recombination  rate  coefficients  are  1.8  x  10  7  cm3/sec  and 
2.0  x  10' /  cm3/ sec  for  and  0^,  respectively.  Therefore,  one  can 
always  compare  this  with  the  attachment  and  diffusion  rates  to  see 
whether  recombination  should  be  included.  The  three-body  attachment 
of  electrons  and  the  dissociative  attachment  with  are  also  tenpera- 
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ture  dependent.  The  three-body  process  has  a  rate  coefficient24  of  1.8 
x  10  cn/Vsec  and  1.0  x  10  ^  cmVsec  when  the  third  bodies  are  mole¬ 
cular  oxygen  and  nitrogen,  respectively.  Such  coefficients  will  give 
a  combined  attachment  rate  at  the  sea  level  of  6  x  10  sec  .  This 

implies  that  the  lifetime  of  a  free  electron  due  to  the  three-body 

-8 

attachment  is  in  1.6  x  10  seconds.  This  should  be  compared  with  a 
lifetime  of  5  x  10  seconds  due  to  dissociative  recombination  for  an 
electron  density  of  ng  =  10^  cm  For  lower  electron  densities, 
lifetimes  due  to  dissociative  recombinations,  became  much  longer.  How¬ 
ever,  when  the  neutral  densities  are  reduced  below  one  atmosphere, 
there  may  be  regions  and  electron  densities  of  interest  where  recombin¬ 
ation  should  be  included.  These  rates  should  also  be  compared  with  the 
diffusion  rate  (see  Equation  3-22)  which  generally  becomes  very  high  at 
very  low  pressures. 
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Fig.  3.3  —  The  ionization  frequency  per  Torr  as  a  function  of  pressure  for  several 
incident  radiation  wavelength  and  fixed  pulse  duration  of  10~6  sec 


3-16 


f0 


380  190 

PRESSURE  [ TORR) 

The  corresponding  effective  breakdown  field  for  the  ionization 
frequencies  of  Fig.  3.3  as  a  function  of  pressure 
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■  Threshold  Breakdown  Power  as  a  function  of  pressure 
radiation  wavelength  and  a  fixed  pulse  duration  of  3  x 


for  several 
10~®  sec 
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Fig.  3.6  —  Threshold  Breakdown  Power  as  a  function  of  pressure  for  several  incident 
radiation  wavelengths  and  a  fixed  pulse  duration  of  3  x  10“ 7  sec  (compare  with  longer 
time  of  Fig.  3.5) 
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